Thermally induced directed currents in hard rod systems 
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. We study the non equilibrium statistical properties of a one dimensional hard-rod fluid undergoing collisions 

04 ' and subject to a spatially non uniform Gaussian heat-bath and periodic potential. The system is able to sustain 

finite currents when the spatially inhomogeneous heat-bath and the periodic potential profile display an appro- 
priate relative phase shift, <j>. By comparison with the collisionless limit, we determine the conditions for the 
most efficient transport among inelastic, elastic and non interacting rods. We show that the situation is complex 
as, depending on shape of the temperature profile, the current of one system may outperform the others. 
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Recently there has been an upsurge of interest in the understanding of non equilibrium systems which even in the absence of 
an applied bias can generate currents. Typical examples are the thermal ratchets and Seebeck ratchet J]]], where an asymmetric 
potential and a non Gaussian noise generate a directed motion. 

Several authors showed that a class of geometrically asymmetric elastic objects undergoing some holonomic constraint and 
coupled to heat baths at different temperatures can rectify thermal fluctuations and thus produce work. The absence of a time- 
reversal symmetry invalidates the standard detailed balance 0,111]. The directed motion of microscopic systems of somehow 
different nature was also studied several years ago by Landauer J(| 13] who considered a bistable potential with an hot heat 
reservoir placed at one side of the potential peak and a cold reservoir on the other side and predicted a directed current of 
particles toward the colder side. This is the so-called the blowtorch effect which has been exploited in Refs. JU, [H] to produce 
, directed currents in periodic non-isothermal systems. 
Q ■ In this paper we study how the presence of interactions, such as excluded volume and inelastic collisions among Brownian 
O | particles, affects the blowtorch mechanism. To the best of our knowledge only the case of overdamped independent particles has 
been analyzed in detail ITolfTlll . 
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The model consists of N impenetrable hard-rods of mass, m, size rj and position x, (f ) (; = 1 , . . . ,A0 evolving on a segment of 
length L according to the dynamics |[T2l IT3ll 

d 2 Xj dxi dV(xi) 



m-^ = -my^ + v^yTfr) &(/) (1) 



We assume cyclic boundary conditions, so that particles crossing with positive velocity the point x = L reenter at the point x = 
\ and viceversa. 

^ . Equation ([TJ is based on the assumption that four kind of forces act on the rods. These are: 

i) the frictional force, —mydxj/dt, proportional to the velocity; 

ii) the gradient of a time-independent spatially periodic potential V(x) = Vocos(2nx/w) of period w, tending to confine the 
particles near its minima; 

iii) the stochastic driving force, mimicking the action of a heath bath with a spatially non uniform temperature profile, has an 
intensity T(x) = T c + 7/,s(x) with s(x) = (1 +tanh[jLisin(27r/wx— <p)])/2 a periodic smooth step-like function between and 
1, which alternates cold T c and warm T c + 7), regions of size w in [Q,L\. The values T c ,Tf, characterize the temperature jump 
amplitude and <p = Itixq/w determines the mutual shift of T(x) and V(x). All temperatures are measured in units such that the 
Boltzmann constant is fes = 1 . As usual, ^,(r) is a zero mean and Gaussian noise with auto-correlation (^,-(f )%j(s)) = 8jj8(t —s); 

iv) finally, the term Y,jfij indicates the resultant of contact impulsive forces acting on i due to the particles j ^ i. namely, the 
rods experience mutual inelastic collisions occurring at contact \xt+\ — x,- = a. After each collision the velocities of a pair 

change according to the rule V; = v,- — (1 + a)(vj — vj)/2 and v'j = vj + (1 + a)(v ( - — V/)/2, where the prime indicates 
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post-collisional values and a is the coefficient of restitution. 



III. THEORY 



In the limit y — > °° and my — ► const, the multiple time scale analysis of ref. lfl4j[l5ll can be extend to the present case to derive 
the follo wing evolution of the one-particle density p(x,t) from the Kramers equation for the phase-space distribution function 

f(x,v,t) M- 



dp(x,t) dJ(x,t) 



where the associated current reads: 



J(x,t) 



dt 



T g [x)p{x,x) -F(x)p(x,t) + 



dx 
(1 + a) 



= 



(2) 



my I dx 

-T g (x-o)g 2 {x,x-o)p(x-o,t) } 



p(x,t) T g {x+a)g 2 {x,x+a)p{x+a,t) 



(3) 



The first term represents the single particle contribution to the current, whereas the second term, non linear in the density, 
accounts for the excluded volume effect. It contains the pair correlation function g2(x,x') evaluated at contact, the coefficient of 
restitution a and the granular (or kinetic) temperature T g of particles related to the heat bath temperature by 



T g (x) = T(x){ 
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gi{x 1 x + a)p {x + a) + g2(x,x - a)p (x - a) j 



(4) 



In the high density limit, the above equations indicate that inelasticity, temperature and density itself result intimately connected 
to determine the system transport properties. On the other hand, in the opposite limit, the latter term in Eq. (0 can be neglected 
and analytical expressions of current J and particle density p (x) can be explicitly worked out. Non interacting system is a 
meaningful comparison as it constitutes the low-density regime extrapolation of interacting particle behavior. The non interacting 
system admits two types of stationary solutions: those corresponding to vanishing and non vanishing current J respectively. The 
key quantity determining the presence of a systematic flux is the "entropy" integral 



S(x) 
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Using the periodicity of the system, we obtain the stationary current Jo 



JS(w) 



J ° Y ac-b[l-^ 



in terms of three constants a, b,c related to S(x) by 



a = dx — — - 
Jo T(x) 



e -s(x) 

dx — — - , 

T(x) Jo 
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S(x) 



(5) 



(6) 



It can be easily verified that, when V(L) = V(0) (zero external load) and T(x) is constant or <j) = 0, the current automatically 
vanishes. It is instructive to discuss how the current of the non-interacting system depends on the temperature scales T c and 7),. 
The current, for T c fixed, grows as the temperature step 7], increases, for the jumps over the barrier to become more prob able. On 
the other hand, if 7), is fixed, J does not depend monotonically on the temperature T c . Clearly, in the limit T c —> 0, the particles 
have a small probability to escape from the potential minima, so that the current must vanish. In the opposite limit T c Vo, the 
confining effect of the energy barriers is negligible and the current must vanish too. For intermediate values of T c , a maximum 
in the current is expected at a r c -value which runs to zero as Vb is reduced. 



IV. NUMERICAL RESULTS 



The case of interacting particles is not so fortunate as not amenable to analytic solution, therefore our study will be based 
mainly on numerical simulations. The picture, indeed, remains qualitatively but not quantitatively similar to that of the non 
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FIG. 1: (Color online) Density profiles of the systems at constant T c /Vq = 0.15, 7/,/Vo = 0.25 and increasing packing fraction tj = 1/8 
(orange), r\ = 1/4 (black), 7] = 1/2 (red) and r] = 3/4 (blue). The dashed curve corresponds to the NI system, while the closed circles and 
open squares symbols correspond to the interacting systems with a = 1 and a = 0.8 respectively. The other parameters are the following: 
m= 1, a = 1, fx = 4.0 and (j) = it. 



interacting case and reveals some interesting features. 

When considering interacting systems, two main effects come into play which may lead to a behavior deviating from the non 
interacting system. First, the mutual repulsion between particles induces dynamical correlations, either promoting the exit from 
a potential well via energetic collisions or forbidding a jump towards a too crowded well. A meaningful parameter commonly 
used to take into account the crowding degree of a granular system is the packing fraction 77 = No /L (0 < 77 < 1). Second, the 
granular temperature of an inelastic system is generally lower than the temperature of the elastic counterpart. Thus, thermally 
activated transport is expected to be less efficient when dissipative collisions are at work. 

One may wonder on the specific influence of the model parameters T c /Vo, T^/Vq, 0, the packing 77 and inelasticity a on the 
transport properties of the inelastic system (Inel). For sake of shortness, here we discuss only temperature effect choosing (j) = K 
(which optimizes the current) and we explore some significant regimes in 77 and a. Moreover we compare the results of the 
inelastic system with the elastic (El, a = 1) and non-interacting (NI, 77 — > a /L) ones, by also analyzing the conditions for the 
efficient transport. 

It is convenient to start the discussion by considering first the stationary density profiles shown in Fig. Q] The inspection of 
profiles, indeed, provides a first indication on the way particles react to parameter variation and how they distribute over the 
effective landscape generated by the potential and temperature profiles. In the NI case, with 7), <C Vrj an d T c <C Vq, the combined 
effect of temperature and potential profiles preferentially confines the particles in the narrow region determined by the minimum 
of the potential and the nearest colder temperature zone; in other words, the phase difference between the minima of T(x) 
and V (x) produces a sort of "cage" trapping the particles with small momenta. This feature is clearly noticeable in the Pni{x) 
structure of Fig. Q] which develops peaks in the cage-region. The large preferential confinement of the NI system is however 
impossible to particles with excluded volume interaction whose density profiles become soon broader and flatter on increasing 
the average packing fraction 77, see Fig.Q] This corresponds to an effective decreasing of the barriers seen by the interacting 
rods and such differences in the density profiles translates into different transport properties. 

We study thus the dependence of the currents JNi^EhJinei on T c (see Fig. |2| at different packings 77 and 7/,/Vb = 0.25 fixed. At 
low values of 77 = 0.125, the current of the interacting systems (green symbols) behaves like the non interacting one (dashed 
black) as expected. Increasing 77, the currents Jei and Jj ne i while remaining very similar to each other, strongly deviate from 
the corresponding Jnj obtained by multiplying for the appropriate number of particles, N, the single particle current J\. The 
currents of the interacting systems start developing a maximum at lower T c (around T c /Vo ~ 0.5 for 77 = 0.25 or T c /Vq ~ 0.2 
for 77 = 0.5). This behavior can be explained by means of the excluded volume effect that reduces the effective height of the 
barriers as the mean packing of the wells increases. As a consequence the maximum of the current vs. T c curve is located at 
smaller temperatures with respect to the corresponding NI case. However, as the packing becomes sufficiently high a second 
effect come into play which changes the above scenario. In this regime, the effective barriers Vo become enough small to be of 
the same order as the temperature step 7), and the particles can escape from the potential well even if T c — > 0. The currents Jei 
and Jj ne i, for 77 = 0.75, show in fact a finite values for small value of T c /Vo (see blue symbols in Fig.[2j\). Moreover, increasing 
T c , the particles become more energetic and the rectifying effect is reduced determining a monotonically decreasing trend of the 
currents. 

Another aspect to consider is the influence of inelasticity on the transport. The curves in Figj2]\ show that, for 77 sufficiently 
high (e.g. 77 = 0.5), the inelastic system (open symbols) becomes more efficient than the elastic one (closed symbols) as long 
as T c /Vq > 0.4. On the other hand, increasing 77, Jj ne [ is larger than Jei for all values of T c (blue symbols). This behavior can 
be explained analyzing the kinetic temperature fields of the two systems. In Fig.[2j}, we show these fields for 77 = 0.5 and two 
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FIG. 2: (Color online) A: particle current versus the rescaled temperature T c /Vo for four different values of the packing fraction T] = 1/8 
(green), T] = 1/4 (black), T] = 1/2 (red) and 77 = 3/4 (blue). Open and full symbols refer to inelastic and elastic system respectively. For 
sake of comparison to the low 17 regime also the NI current is plotted (dashed curve). B: granular temperature profiles of the systems for two 
different value of T c /Vq: 0.7 (top) and 0.15 (bottom). Open and full symbols correspond to inelastic and elastic system respectively. The 
dashed curve is the temperature profile for a NI system. 7/,/Vrj = 0.25 and the other parameters are the same as in Fig. [T] 



different values of T c /Vq. 

While the elastic and non-interacting profiles of T(x) are very close, the effective value of the temperature, in the inelastic 
case, decreases near the potential minima. This indicates that the ratio T c /Vq is smaller if a < 1. The inelastic system can thus 
be approximately considered as the elastic one, but with a lower temperature T c , i.e. Ji, le i(T c ) « Jei{T c — 8T C ) for a fixed T c . 
Such decreasing of T c helps to rectifying the fluctuations when the energy is enough for the activated thermal process, i.e. for 
77 = 0.75 or T7 = 0.5 and T c /Vq > 0.4, or it reduces the transport when T c + T h <C % (i.e. for T7 = 0.5 and T c /Vq < 0.4). 



V. CONCLUDING REMARKS 



In this work, we have shown that contrary to intuition, one-dimensional systems with hard-core interactions can display a 
more efficient "blowtorch" effect than a non interacting one. Our simulations show this efficiency inversion to be correlated 
to the resistance of hard-rods to localize in narrow regions. Under some conditions, the inelasticity makes the transport more 
efficient by reducing the average kinetic energy. 
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